Evidence for Radio Detection of Extensive Air Showers 
Induced by Ultra High Energy Cosmic Rays 
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Firm evidence for a radio emission counterpart of cosmic ray air showers is presented. By the 
use of an antenna array set up in coincidence with ground particle detectors, we find a collection 
of events for which both time and arrival direction coincidences between particle and radio signals 
are observed. The counting rate corresponds to shower energies > 5 x 10 16 eV. These results 
open overwhelming perspectives to complete existing detection methods for the observation of ultra 
high-energy cosmic rays. 
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The origin and the nature of ultra high-energy cosmic 
rays observed above 10 19 eV is one of the most chal- 
lenging questions in astroparticle physics 0. It is also 
of major importance for the understanding of extremely 
energetic phenomena in the Universe. In the 1960's, ra- 
dio emission associated with the development of exten- 
sive air showers was predicted and measured 0, 0] . A 
flurry of experimental investigations, carried out at that 
time with a wide variety of setups and in different fre- 
quency domains, have provided initial information about 
signals from 10 17 eV cosmic rays 4J. Plagued by difficul- 
ties (poor reproductibility, atmospheric effects, technical 
limitations) efforts ceased in the late 1970's to the bene- 
fit of ground particle [l| and fluorescence detections 
With the advent of cosmic ray research involving giant 
surface detectors as in the Auger experiment Q, the ra- 
dio detection, having potentially 100 % duty cycle and 
sensitive to the longitudinal development of the showers, 
is now reconsidered. In recent years, making use of up to 
date techniques in transient analysis, several groups have 
undertaken the task of reinvestigating the phenomenol- 
ogy of radio pulses H, |j| which are to a large extent 
a terra incognita. Using our experiment CODALEMA 
(COsmic ray Detection Array with Logarithmic Electro- 
Magnetic Antennas) , located at the Nancay Radio Obser- 
vatory, firm evidence for a radio emission counterpart of 
cosmic ray air showers is presented. From 60 days of run- 
ning time, 18 events show both time and arrival-direction 
coincidences between particle and radio signals. This de- 
tection rate corresponds to shower energies > 5 x 10 16 eV. 

During a first phase [ll|, begun in March 2003, the 
trigger was generated using signals from one antenna de- 
voted to this purpose. A good knowledge of the transient 
radio environment was acquired |12J. Coincident tran- 



FIG. 1: CODALEMA setup for the second phase. The par- 
ticle detector acts as a trigger with a fourfold coincidence 
requirement. 



sients were measured from antennas spaced by several 
tens of meters, and information on the arrival direction 
and on the field amplitudes was obtained 01 • A first at- 
tempt to identify air shower candidate events was made 
based on gross features associated with air showers in the 
10 17 eV range [b|. However, due to the rather limited 
knowledge available concerning radio transients from the 
sky in the nanosecond range, the results obtained, al- 
though encouraging, could not be conclusive. 

For the second phase [l^, to make sure that the mea- 
sured transients were associated with air showers, four 
stations of particle scintillator detectors acting as a trig- 
ger have been added. This setup, including a new an- 
tenna configuration, is shown in Fig. ^ It uses seven log- 
periodic antennas of the t ype constituting the Nancay 
DecAMetric array (DAM) [lq. The antenna waveform 
signals are recorded after RF signal amplification (1-200 
MHz, gain 35 dB) by LeCroy digital oscilloscopes (8- 
bit ADC, 500 MHz sampling frequency, 10 /xs recording 
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FIG. 2: Example of a seven antenna event, numerically fil- 
tered in the frequency band 37-70 MHz, obtained in coinci- 
dence with a particle trigger (first trace). 



FIG. 3: Time delay between the radio plane front and the 
particle plane front 



time). To get enough sensitivity with these ADCs, the 
antennas are band pass filtered (24-82 MHz). 

The trigger corresponds to a fourfold coincidence 
within 600 ns from the particle detectors originally de- 
signed as prototype detectors for the Auger array [t7| . 
Each 2.3 m 2 module (station) has two layers of acrylic 
scintillator, read out by a photomultiplier placed at the 
centre of each sheet. The photomultipliers have copper 
housings and it has been carefully verified that no corre- 
lation exists between high photomultiplier signal ampli- 
tudes and the presence of antenna signals. The signals 
from the upper layers of the four stations are digitized (8- 
bit ADC, 100 MHz sampling frequency, 10 /xs recording 
time). A station produces a signal when a coincidence 
between the two layers is obtained within a 60 ns time 
interval. This results in a counting rate of around 200 Hz. 
The rate of the fourfold coincidence of the four stations 
is about 0.7 event per minute. The stations are located 
close to the corners of the DAM array, and in this con- 
figuration the particle detectors delimit an active area of 
roughly 7 x 10 3 m 2 . Using arrival times from the digi- 
tized photomultiplier signals, it is possible to determine 
the direction of the shower by triangulation with a plane 
fit. From the arrival direction distribution, a value of 
16 x 10 3 m 2 xsr is obtained for the acceptance, which 
corresponds to an energy threshold of about 1 x 10 15 eV. 

For each fourfold coincidence from the particle detec- 
tors, the seven antenna signals are recorded. Due to the 
relatively low energy threshold, only a small fraction of 
these air shower events is expected to be accompanied by 
significant radio signals. 

The recognition of the radio transients is made during 
the offline analysis (see Fig. EJ. For each waveform, the 
maximum value within the 37-70 MHz band is searched 



for in the time interval [2 /js ; /xs]. Then the noise 
level is determined using the mean /x„ and standard de- 
viation a n of the signal outside of a ±1.2 /xs zone around 
this maximum. A pulse transient is flagged as "detected" 
when, in a ±0.3 /xs zone around the maximum, the av- 
erage of the squared signal /i s satisfies the condition 
Ms > l^n + k x a n where the number k is adjusted empir- 
ically to optimize the rejection of non pulse-like signals. 

For each antenna, if the threshold condition is fulfilled, 
the time of the maximum is determined. When at least 
3 antennas are flagged it becomes possible to apply a 
triangulation procedure |14| , and the event is declared to 
be a radio candidate if the arrival direction obtained is 
plausible, i.e. above the horizon. 

On average, there is one radio event for every 10 3 par- 
ticle triggers. This indicates that the energy threshold 
for radio detection is substantially higher than that of 
the particle array. 

Event selection then proceeds using the arrival times. 
The time of passage of the radio wavefront through a 
reference point is compared to the particle front time ex- 
tracted from the scintillator signals. The distribution ob- 
tained is shown in Fig.|3J The data correspond to 59.9 ac- 
quisition days and 70 antenna events. A very sharp peak 
(a few tens of nanoseconds) is obtained showing an unam- 
biguous correlation between certain radio events and the 
particle triggers. This peak is centred around 40 ns with 
a systematic error of 20 ns. The delay between electric 
field and particle shower times could be measured with 
our apparatus, though more thorough studies and higher 
statistics will be necessary. As expected, in the 2 /is 
window where the search is conducted, radio transients 
can occur which are not associated with air showers and 
correspond to accidentals. Uncorrelated to the particles, 
such events have an uniform arrival time distribution. 
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FIG. 4: Event arrival directions. The circles indicate direc- 
tions reconstructed from antenna signals whereas the crosses 
correspond to directions given by the scintillators. Except for 
the marked event, each circle is associated with the nearest 
cross. 

Finally, if the time-correlated events correspond to ex- 
tensive air showers, the arrival directions reconstructed 
from both scintillator and antenna data should be 
strongly correlated. Fig. 0] exhibits, for the 19 events 
located in the main peak of Fig. |3 the reconstructed ar- 
rival directions. The distribution of the angle between 
the two reconstructed directions is, as expected, a gaus- 
sian distribution centred on zero multiplied by a sin func- 
tion coming from the solid angle factor. The width of the 
corresponding gaussian is about 3.5 degrees. The angle 
for the remaining uncorrelated event is much bigger, thus 
this event is certainly an accidental: its arrival direction 
given by the antennas is close to the horizon, which is 
typical of events from radio interference due to human 
activity. Moreover, the presence of one chance event in 
the peak is fully compatible with the observed uniform 
distribution in the 2 /is window. 



For the 18 correlated events of Fig.0]the antenna mul- 
tiplicity varies from 3 to 7 and the amplitude of the fil- 
tered signals goes up to 1.2 mV/m. In the case of the 
seven-antenna event of Fig. [3 all the filtered signal am- 
plitudes are around 0.25 mV/m. This signal level corre- 
sponds to typical values expected for air shower energies 
in the 10 17 eV range In order to obtain a rough 
estimation of the energy threshold for our setup, we sup- 
pose that the acceptance is the same for the two types 
of detectors (i.e. 16 x 10 3 m 2 xsr). The observed event 
rate then leads to an approximate energy threshold of 
5 x 10 16 eV. 

In conclusion, these results strongly support the claim 
that electric field transients are generated by extensive 
air showers and that they can be clearly measured us- 
ing wide band techniques. More data and technical up- 
grading are needed to examine the contribution that ra- 
dio detection could bring in determining the energy and 
the nature of the cosmic rays. Various improvements 
are planned for CODALEMA: setting up of additional 
scintillators will make possible the determination of the 
shower energy and core position; extension of the W-E 
antenna line will enable to better sample the radio sig- 
nal spread; increase of the ADC dynamics up to 12-bit 
encoding will allow to record the full 1 - 100 MHz fre- 
quency band, shower parameters could then be inferred 
from the full signal shape. In a subsequent upgrade, it 
is planned to install autonomous dipoles equipped with 
active front-end electronics, self-triggered and self-time- 
tagged. This is part of current investigations to use the 
potentialities of the radio detection techniques in con- 
junction with existing surface detectors. Among various 
particles capable to generate air showers, neutrinos are 
thought to be probably the only ones which could address 
the hadronic nature and source location of ultra energetic 
cosmic rays. Taking into account the ability of the ra- 
dio detection method to characterize inclined extensive 
air showers 0] , the presented results could considerably 
enlarge the scope of such studies for cosmology related 
phenomena. 
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